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SUMMARY 

The  unsteady  transonic  small  perturbation  differential  equation  is  converted  into  an 
integro-differential  equation  by  tq/plication  of  the  classical  Green’s  function  method.  It  is 
shown  that  no  contribution  from  shock  waves  explicitly  appears  in  this  integral  equation, 
due  to  the  shock  capturing  properties  of  the  Green's  function  method.  After  assuming  that 
the  motion  consists  of  small  infinitesimal  perturbations  around  a  thin  nearly-planar  body, 
a  simplified  integral  equation  for  the  streamwise  velocity  component  is  obtained,  which  is 
suitable  for  fast  numerical  computations. 
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pressure  coefficient 

scaled  pressure  coefficient  for  A/  <  1 

scaled  prmure  coefficient  for  M  >  1  _ 

domain  function 

Green’s  function 

Heaviside  step  function 

coefficients  for  subsonic  case 

coefficients  for  supersonic  case 

typical  chord  length 

free  stream  Mach  number 

unit  normal  in  Prandtl-Glauert  variables 

unit  normal  in  Eulerian  variables 

{(x-A:i)2-«*(0>-yi)2+(z-zi)*l}» 

xi)2 +/S2[0'->-i)2 +(z-zi)*]}* 

{(X-Xi)^-i  Y~  Ki)2-(Z-Zi)*}‘ 
projection  of  wing  surface  onto  z  =  0  plane 
projection  of  wake  surface  onto  z  =  0  plane 
equation  defining  wing,  wake  and  shock  surfaces 
stream  wise  perturbation  velocity  {=d<bldx) 

=  dOp/dx 

free  stream  velocity 
Eulerian  variables 
Prandtl-Glauert  variables 
position  of  trailing  edge 

ratio  of  specific  heats 
2-(2-y)A/* 

thickness  ratio,  Dirac  delta  function 
acoustic  time  delay  between  two  points 
acoustic  time  delay  in  Prandtl-Glauert  variables 


\x-Xi\ -{(y-  yi)*+(z-Zi)»}* 
free  stream  density 

surface  surrormding  wing,  wake  and  shock  waves 
perturbation  potential 
scaled  perturbation  potential 
^solution  for  purelyLsymiDetrical.  flow  . 

velocity  potential 


1.  ii^ODucno^  .  _ _  _  ^  _  _  .  _ 

The  integral  equation  method  in  steady  transonic  flow  problems  was  first  used  by  Oswatitsch 
[13].  Since  then,  the  method  has  been  extended  and  modified  by  Spreiter  and  Alksne  [IS], 
Noistrud  [11],  Nixon  [9]  and  Piers  and  Sloof  [14].  Nixon  [8],  [10]  extended  the  integral  equation 
method  to  harmonically  oscillating  two-dimensional  airfoils,  by  perturbing  from  the  steady  solu¬ 
tion.  A  weakness  of  this  formulation  is  its  inability  either  to  create  or  eliminate  shock  waves  as 
a  result  of  the  unsteady  motion.  In  this  paper  we  use  a  method  initially  developed  by  Morino  [6], 
[7]  for  purely  subsonic  and  supersonic  oscillatory  flows.  The  method  uses  the  classical  Green’s 
theorem  approach  to  derive  an  integro-difierential  equation  for  the  perturbation  velocity  potential 
for  unsteady  flow  about  a  general  body.  Recently  the  technique  was  adapted  to  unsteady  tran¬ 
sonic  flow  by  Tseng  and  Morino  [16]. 

In  the  following,  the  isentropic  invisckt  flow  of  a  perfect  gas,  initially  iirotational,  is  con¬ 
sidered.  Under  this  hypothesis,  the  existence  of  a  velocity  potential  may  be  assumed  such  that 
the  fluid  velocity  v  =  (u,  v,  w)  =  v^.  The  most  basic  approximation  of  inviscid  aerodynamics 
is  that  of  a  small  disturbance.  Thus  if  £/<»  is  the  free  stream  velocity  (which  is  assumed  to  be 
directed  solely  in  the  x-direction),  then  a  perturbation  potential  ^  can  be  defined  such  that 

y,  z\  I)  =  Un[x+4(x,  y,  z;  /)],  (1) 

where  x,  y,  z  represent  a  rectangular  cartesian  co-ordinate  system  and  /  is  the  time  variable. 
Note  that  subsequent  equations  are  expressed  in  non-dimensional  co-ordinates  based  on  a  length 
scale  /,  a  typical  value  of  the  airfoil  chord  length,  a  velocity  scale  Ua>,  and  a  density  scale  pm, 

a  typical  value  of  the  density  at  infinity.  The  time  variable  is  then  scaled  with  l/Um  and  the 

pressure  with  pm  Um^. 

In  terms  of  the  dimensionless  Eulerian  co-ordinates  x,  y,  z  and  t,  the  governing  equation 
for  the  perturbation  velocity  potential  can  be  written  as 

Id  a\* 

where  M  is  the  free  stream  Mach  number  and  %  includes  all  the  nonlinear  terms.  For  simplicity 
the  discussion  will  be  restricted  to  the  nonlinear  term 

X  =  A/*(l+y*)^x^x*,  (3a) 

with 

y*  =  2-(2-y)A/2,  (3d) 

where  y  is  the  ratio  of  specific  heats.  The  procedure  used  to  derive  the  approximate  equation 
[(2)  with  X  given  by  (3fl)],  when  the  flow  is  transonic,  is  based  on  the  small  parameter  8  repre¬ 
senting  the  ratio  of  airfoil  thickness  to  chord  length.  The  main  assumptions  are  that  as  8  -»  0, 
the  quantities  (1  —  A/*)/«(8),  ylp{S),  zlp{d)  and  f//n^8)  all  remain  fixed,  whereas  «(8), 

/i(8)  0  (cf.  with  [I]  and  [2]).  If  terms  smaller  than  tp*  and  e®  are  neglected  then  y  is  given  by 

equation  (3a).  Note  that  the  coefficient  of  ^  ^xx  in  (3a)  is  usually  written  as  I  +  y).  The 
modification  in  (3a)  comes  forward  if  the  small  perturbation  equation  is  derived  directly  from 
the  conservation  of  mass  equation  (see  [17]). 

The  effect  of  this  modification  is  that  the  relations  for  a  normal  shock  show  better  agree¬ 
ment  with  those  of  full  potential  theory.  It  also  leads  to  a  better  agreement  of  the  critical  pressure 
coefficient  with  its  exact  isentropic  value  (i.e.  when  the  steady  local  fluid  speed  equals  the  local 
sonic  speed)  (see  [3]  and  [4]). 


I 


In  the  subsequent  analysis  a  general  theory  of  potential  aerodynamic  flow  around  thin 
lifting  bodies  having  arbitrary  shape  and  motion  is  presented.  In  Section  2,  for  subsonic  free 
stream  Mach  numbers,  the  small  perturbation  differential  equation  (2),  is  converted  into  an 
integro-difierential  equation  for  the  velocity  potential,  by  application  of  the  classical  Green’s 
theorem  approach.  It  is  also  shown  that  no  contribution  from  shocks  explicitly  appears  in  the 
subsequent  integral  equation  for  due  to  the  shock  capturing  nature  of  the  method.  In  Section  3 
and  4  the  intejp^  equation  formulation  is  simplified  by  assuming  t^at  the  motion  consists  of 
sii^riniihit^mal' perturbations  around  a  stcEufy  thin  n^rly-planar  body.  Eventually  we  obtain," 
a  simplified  integral  equation  for  the  streamwise  perturbation  velocity  component,  which  is 
suitable  for  fast  numerical  calculations.  In  the  Appendix  a  similar  integral  equation  is  derived 
for  the  case  when  the  free  stream  Mach  number  is  supersonic. 


2.  SOLUTION  BY  NONLINEAR  GREEN’S  FUNCTION  METHOD 

The  method  of  solution  for  (2)  is  based  upon  the  well-known  Green’s  function  technique. 
The  Green’s  function  for  the  equation  of  potential  is  the  solution  of  the  problem 


/a  av 

Va  ^  =  Kx-xi,y-yu  z-zu  t-h) 


with 


G  =  0  at  infinity. 


(4o) 


(4b) 


where  b  is  the  Dirac  delta  function.  The  solution  of  (4a  and  4b)  for  the  subsonic  case  (3f  <  1) 
is  given  by  (see  Morino  [7]) 


C(x-xi,  t—h)  = 


I 


4»rr^x,  xi) 


where 


and 


r^x,  xi)  =  {(x-xi)a+/3a[(y-;>i)a+(z-2i)a]}», 

jS2=  I -4/a 


(5a) 

(5b) 

(5c) 


®  xi)-A/(x-xi)]. 


(5d) 


The  corresponding  solution  for  the  supersonic  case  (W  >  I)  is  given  in  the  Appendix. 

Using  the  procedure  developed  in  [7]  in  order  to  obtain  a  representation  of  the  perturbation 
potential  in  terms  of  its  value  and  the  values  of  its  derivatives  on  the  wing  body,  the  wake  and 
the  shock  waves  (if  any)  we  find  it  convenient  to  define  a  domain  function  E  such  that 

outside  E 

on  E  (6) 

inside  E, 


where  E  is  a  surface  which  completely  surrounds  the  wing,  its  wake  and  any  shocks,  and  which 
is  defined  by  the  equation 


Following  the  general  Green’s  theorem  method  we  now  multiply  the  equation  of 
aerodynamic  potential  (equation  (2))  by  the  Green’s  function  G  and  subtract  equation  (4a) 
multiplied  by  After  suitable  manipulations  we  find  that 


d/  M  .dG\ 

. 


=  Gx-<l>S(x-xuy-yi,  z-zi;  t-h). 


(8a) 


with 


d  d  ^ 
dt  8t^dx 


(8b) 


In  (8a)  the  arguments  of  <t>  and  its  derivatives  are  xi,  yi,  zi  and  h,  while  the  arguments  of  G  and 
its  derivatives  are  x—xi,  y—yi,  z—zi  and  t—ti.  If  equation  (8a)  is  multiplied  by  the  domain 
function  E,  and  integrated  over  the  whole  four-dimensional  space-time,  it  is  found  that  the  poten¬ 
tial  satisfies  (see  Morino  [6],  [7]), 


00 


4irE(x;  f)^x; 


—  00 


111 


Ifl 


d£»(xi) 


j  I  [v., . 


(9) 


In  equation  (9)  the  superscripted  0  indicates  evaluation  at  time  fi  =  l—O,  where  0  is  defined  by 
(Sd),  thus  S*  indicates  the  surface  defined  by  S®  =  S(jri, >’i,  zii  1—0)  =  0.  Therefore  2*  is  a 
surface  of  the  three-dimensional  space  (xi,  yu  zi),  which  depends  parametrically  upon  x,  y,  t 
and  t.  Note  that  the  surface  2  is  composed  of  three  branches.  The  first,  Eb,  is  the  surface  of 
the  body.  The  second  is  the  surface  2w,  of  the  wake.  Note  that  this  surface  Ew  is  considered 
twice,  since  E  is  a  closed  surface.  In  other  words,  the  two  sides  of  the  wake  are  considered  to  be 
independent  surfaces  having  the  same  equation  but  opposite  outwardly  directed  normals.  The 
third  one  is  the  surface,  Es,  of  any  shock  waves,  for  which  similar  considerations  are  valid. 
In  the  subsequent  analysis  it  is  convenient  to  isolate  the  contribution  of  each  of  these  surfaces. 

If  we  now  substitute  for  x  in  (9)  with  its  approximate  value  (3a),  we  see  that  the  volume 
integral  on  the  right  of  (9)  can  be  integrated  by  parts  in  the  x-direction.  This  integration  by 


parts  delivers  integrals  over  the  bounding  surfaces  of  the  volume.  If  we  separate  the  contribution 
of  each  of  these  surfaces  it  is  found  that  (9)  reduces  to 


00 


47rEix,  t)^x-,  t)  = 


M\1  +  y* 
-  2 


—  00 


(Sb  +  Sw)® 


(Sb+£w)* 

(2b  +  2w)® 

The  contribution  from  shocks  does  not  explicitly  appear  in  (10),  this  is  due  to  the  shock  boundary 
conditions. 


A^  =  0, 


(llfl) 


-A/2A^^x+^t)^+ a|^^i  -Af2^x- 

where  A /  represents  the  jump  in  / across  the  shock.  The  condition  (1  lo)  follows  from  the  fact 
that  there  can  be  no  circulation  around  infinitesimal  paths  threading  the  shock  front,  while  (1  Ih) 
is  found  by  integrating  the  conservation  or  divergence  form  of  (2)  (with  x  given  by  (3a))  across 
a  jump  discontinuity.  The  lack  of  an  integral  over  the  surface  Sg,  does  not  necessarily  signify 
that  shocks  do  not  contribute  to  the  potential  <f>.  As  shown  by  Tseng  and  Morino  [16]  the  con¬ 
tribution  from  shock  waves  is  now  effectively  imbedded  in  or  captured  by  the  volume  integral 
in  (10). 


3.  QUASI-STEADY  APPROXIMATION 

In  practical  applications  it  is  convenient  to  introduce  a  restrictive  assumption,  which  allows 
a  considerable  simplification  of  equation  (10).  In  particular  the  unsteady  motion  of  the  aircraft 
can  be  assumed  to  consist  of  small  infinitesimal  perturbances  around  the  steady  state  configura¬ 
tion,  so  that,  for  the  purpose  of  evaluating  the  integrals  in  (10),  £  and  S  are  assumed  to  be  fixed. 


However,  when  the  boundary  conditions  on  the  body  are  applied,  the  time  dependence  of  the 
surface  position  must  be  taken  into  consideration,  in  order  to  produce  any  unsteady  motion. 
It  is  also  convenient  to  introduce  the  Prandtl>Galuert  variables 


I 


,1 


I 


and 


with 


^  =  jc,  Y  =  py,  Z  =  fiz 


02a) 


(126) 


Substituting  (12o  and  126)  into  (10)  and  applying  the  quasi-steady  approximation,  it  is  found 
that  the  perturbation  potential  satisfies 


£(X)d»(X,  T) 


8ir 


I  S;  k  ““Ik . 


—  00 


-iJJ 


2b+Sw 


"Hr 

£b  +  2w 


(S-i) 

T^|X-Xl|' 


;«fS(Xl) 


I 

'*'4^ 


J — !?li  rfz{Xi) 

I  8/VVlX-Xil/  IX-Xiia/VjlTi-T-e  ^ 


(13a) 


2^b+2Iw 


where 


©  =  IX-Xil-A/lJr-A'i), 

|X-Xi!  =  {(A'-A'i)*-l-(  Y-  Ki)2  +  (Z-Zi)*}«, 

t/  =  'I>x 


(136) 

(13c) 

(13rf) 


and 


dN 


/=  ^/. 


Here  4  and  Si  are  the  gradient  operator  and  unit  normal  to  the  surface  in  PrandtI-Galuert  space. 
It  should  be  noted  that  all  the  terms  on  the  right  hand  side  of  (13a)  have  a  physical  interpretotion 
analogous  to  potential  theory.  The  first  integral  on  the  right  (of  (13a),  which  originates  from  the 
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nonlinearity  on  the  right  of  (2),  has  the  form  of  a  distribution  of  streamwise  directed  doublets 
of  strength  (IP(^  T)/9it  per  unit  volume.  The  second  and  third  integrals  represent  distributions 
of  soiuces  over  the  wing  and  the  wake  with  strengths  given  by  —Xl^nd^jdN  and  (2A/Ot+  t/®/2) 
(f) .  i)/4ir  per  unit  area.  The  last  integral  on  the  right  of  (13a)  corresponds  to  a  distribution  of 
doublets,  directed  normally  to  the  wing  and  the  wake,  with  strength  <h(X,  r)/4ir  per  unit  area. 
Hence  (13a)  asserts  that  the  disturbance  outside  the  closed  surfaces  £b  and  Sw>  is  the  same 
disturbance  that  would  be  produced  by  a  fictitious  distribution  of  sources  and  doublets.  Also, 
it  should  be  noted  that  0^  (see  (13b))  Is  equal  tb  the  time  necessary  for  a  disturbance  to  propogate 
from  Xi  to  X  at  the  speed  of  sound.  In  other  words,  0  is  equal  to  the  usual  acoustic  time  delay. 

At  this  stage  in  the  analysis  we  should  note  that  the  volume  and  surface  integrals  in  (1 3a) 
have  singular  contributions  when  Xi  equals  X.  In  the  preceding  and  subsequent  analysis  it  is 
assumed  that  these  singular  points  are  excluded  from  the  integration  region  by  suitable  principal 
value  definitions.  For  the  volume  integral  we  define  a  principal  value  by  surrounding  the  singu¬ 
larity  with  a  sphere  of  small  radius  and  take  the  limit  as  the  radius  tends  to  zero.  When  X  is  on 
either  the  body  or  the  wake  the  i  he  surface  integrals  are  singular  when  Xi  equals  X.  In  this 
case  we  define  the  principal  v  .iue  by  surrounding  the  singularity  with  a  circle  of  small  radius 
and  take  the  limit  as  the  radtut;  tends  to  zero. 

The  integral  equation  for  the  perturbation  potential  can  be  further  simplified  by  applying 
the  relevant  boundary  conditions  on  the  body  and  the  wake.  Firstly  we  consider  the  kinematic 
condition  on  a  surface  in  an  invisied  fluid.  This  condition  amounts  to  the  statement  that  a  fluid 
particle  initially  in  a  surface,  remains  in  that  surface  throughout  the  motion.  If  the  surface  is 
defined  by  s(x;  t)  —  0  then  the  kinematic  condition  specifies  that 

D 

— (S(x;  /))  =  St  +  Sx  +  ^x-Sx  +  ^y-Sy  +  ^ziSz  =  0 


on 


Six;  t)  =  0, 


(14) 


where  D/Dt  is  the  Eulerian  derivative.  Now  consider  df^ldN,  the  normal  derivative  of  ‘h  in 
Prandtl-Glaucrt  variables.  Using  (12a  and  1 2b),  d<^ldN  can  be  written  as 


bO  ^  MH\+Y*)d<f>,  |VS| 

V- 1  ''rv5i 


(15) 


where  b^/bn  is  the  normal  derivative  of  <t>  in  Eulerian  coordinates,  and  is  known  from  (14)  when 
the  surface  S  is  specified. 

In  the  wake  the  usual  Kutta  condition  applies.  The  wake  is  a  surface  of  discontinuity  in 
where  the  two  sides  of  the  wake  are  independent  surfaces  having  the  same  equation.  If  the  boun¬ 
dary  condition  (14)  is  applied  to  both  sides  of  the  wake  we  have  that 


on 


Sw(x,/)  =  0,  (16) 

where  A(  )  now  represents  the  value  on  the  upper  side  minus  the  value  on  the  lower  side  and 
5w  specifies  the  wake  surface.  Equation  (16)  simply  states  that  there  is  no  jump  in  normal 
velocity  across  the  wake.  Also  from  the  Kutta  condition  we  know  that  the  flow  leaves  the  trailing 
edge  smoothly.  Consequently  there  can  be  no  discontinuity  in  pressure  at  the  trailing  edge  and 
throughout  the  wake.  When  the  flow  is  isentropic  the  appropriate  approximate  expression  for 
the  pressure  coefficient  is 


Cp  =  — 2(^x+^t). 


(17) 
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Thus  across  the  wake  we  have  that 


A(^x) + Mfh)  —  0 


on 


Svfix,  r)  =  0 
or  in  Prandtl-Glauert  variables 
A/2 

-  — Ai/  = 


(18) 


on 


Sw(X,  T)  =  0. 


(19) 


If  we  now  apply  (16)  and  (19)  to  (15)  we  find  that  across  the  wake, 

a(^)  =  (51 .  i)A/A<l)T 

on 

Sw(X,  T)  =  0.  (20) 

Equation  (13a)  can  now  be  simplified  by  substituting  (IS)  and  (20)  into  the  second  integral 
on  the  right  hand  side  (i.e.  (15)  is  used  over  Sb  and  (20)  over  Sw).  However,  it  is  convenient 
at  this  stage,  due  to  the  discontinuity  in  <t>,  to  assume  that  the  wake  can  be  represented  by  a  vortex 
sheet  parallel  to  the  A'T-plane.  Thus  the  normal  direction  to  the  wake  is  parellal  to  the  Z-axis 
and  terms  involving  fl .  i  (in  (13a)),  one  zero  over  the  wake.  Also,  since  the  fluid  is  all  moving  at 
approximately  the  same  stream  speed,  a  vortex  element  at  a  general  point  X  of  the  wake  was 
shed  in  the  past  at  a  moment  determined  by  the  time  interval  T))/A/  required  for  it  to 

reach  X. 

Hence  at  a  general  wake  point  X  and  time  T. 


Af(X,  Y,  Z,  T)  =  Af(XT{  K),  Y,  Z,  T-^{X-  Xt{  Y))IM),  (21) 


where  A't(  Y)  specifies  the  position  of  the  trailing  edge  and /in  (21)  represents  <I>  and  its  derivatives. 
Thus  using  (21)  integrals  over  the  wake  can  be  evaluated  once  <1>  is  known  at  the  trailing  edge. 

Now  using  (15)  and  (20),  in  the  second  integral  on  the  right  of  (13a),  together  with  (21)  on 
the  horizontal  wake,  we  find,  after  a  few  minor  manipulations,  that  the  integro-difierential 
equation  for  O  becomes 


£(X)0(X,  T)  =  <t>p(X,  T) 
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where 


and 


Op(X,  D  =  - 


A/2(l  +  y*) 


^(xi.  t- 
dn 


■«) 


IViSI 


IX-Xil  IviSI 


^,'/2:(Xi) 


Cp  =  -‘HU+p^jM)  =  MHl  +  Y*)CplP^ 
r  =  0+/32(^i-A'T<>'i))/A/. 


(22a) 


(22*) 

(22c) 

(22d) 


In  (22a)  Op  represents  the  solution  for  the  purely  symmetrical  problem  of  linearized  flow  past  a 
thin  symmetrical  wing  at  zero  incidence.  Once  the  shape  of  the  body  is  specified  Op  is  completely 
determined  by  (22*)  and  (14).  In  the  last  integral  of  (22a)  R\t  represents  the  projection  of  the 
wake  surface  Sw  onto  the  X—Y  plane,  which  is  assumed  to  be  the  mean  horizontal  surface  of 
the  wake.  Also,  we  note  here  that  F,  see  (22d),  is  the  time  necessary  for  a  disturbance  to  be  ad- 
vected  downstream  from  the  trailing  edge  to  the  wake  point  Xi  at  free  stream  speed,  plus  the 
acoustic  time  delay  between  Xi  and  the  field  point  X. 

In  order  to  solve  for  the  potential  it  is  necessary,  in  general,  to  obtain  a  numerical  approxi¬ 
mation  for  (22n).  We  note  that  (22a)  determines  the  potential,  at  time  T,  on  the  body,  the  wake 
and  in  the  fluid  from  known  values  of  <*,  <t>T  and  V  at  times  less  than  T.  The  evaluation  of  <I>t 
and  U  from  <I>  involves  the  use  of  finite  differences.  When  evaluating  U,  due  care  must  be  used  to 
ensure  that  the  correct  differencing  scheme  is  used  in  regions  where  the  flow  is  locally  super¬ 
sonic.  It  should  be  noted  however,  that,  the  necessity  to  find  U  by  finite  differencing  can  be 
eliminated  by  differentiating  (22a)  with  respect  to  X  (taking  due  care  of  the  singularity  in  the 
volume  integral).  Thus  equation  (22n)  and  the  integro-differential  equation  for  U  can  then  be 
used  with  time  differencing,  to  step  forward  in  time  evaluating  <!>  and  V  at  all  points  of  the  flow. 
However  at  this  stage  it  is  worthwhile  making  a  further  approximation  which  in  effect  assumes 
that  the  wing  is  a  thin  nearly  planar  body. 


4.  THIN  BODY  APPROXIMATION 

The  determination  of  the  original  equation  (2)  (with  (3a))  was  based  upon  the  smallness 
of  the  parameter  S  (which  represents  the  ratio  of  airfoil  thickness  to  chord  length).  Consequently, 
it  is  consistent  to  assume  that  the  wing  is  a  thin  nearly-planar  body  with  normal  direction  parallel 
to  the  Z-axis.  Applying  this  approximation  to  the  surface  integrals  in  (22a),  but  not  Op,  produces 
a  significant  simplification.  The  third  term  on  the  right  of  (22a)  becomes  zero  and  the  last  two 
integrals  may  be  integrated  by  parts  in  the  Z-direction,  to  produce  an  integral  over  a  vortex 
sheet.  Note  that  this  is  the  same  as  analytically  continuing  the  integrands  through  the  surface 
of  the  body  onto  a  mean  horizontal  plane  and  then  ignoring  terms  of  second  order  in  *. 


Applying  the  thin  body  approximation  to  the  surface  integrals  in  (22a)  and  subsequently 
integrating  by  parts  the  last  two  integrals,  we  find  that  the  integro-differential  equation  for  can 
be  written  as 


^(XWX,  T)  =  Op(X,  7^ 
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U(XTiYi),Yifi,T-r), 


(23) 


where  Rb  represents  the  projection  of  £b  onto  the  X~  Y  plane  and  (19)  has  been  used  to  eliminate 
A<hT  on  the  wake.  Note  that  the  last  two  integrals  in  (23)  represent  lifting  effects,  and  correspond 
to  a  vortex  distribution  over  the  A'—  y  plane.  It  should  also  be  noted  that  if  the  flow  is  steady, 
then  derivatives  with  respect  to  T  are  zero,  A{/  is  zero  at  the  trailing  edge  and  throughout  the  wake, 
and  (23)  then  reduces  to  the  formula  derived  by  Klunker  [5]. 

The  integral  equation  (23)  for  the  reduced  velocity  potential  contains  derivatives  of  the  re¬ 
duced  velocity  potential  under  the  integral  signs  and  as  such,  is  an  integro-differential  equation 
for  <I>.  Although  it  is  quite  conceivable  to  solve  this  equation  by  straight  forward  numerical 
means  it  is  convenient  to  derive  an  integral  relationship  for  the  streamwise  velocity  component  U. 
The  relevant  integral  equation  is  found  by  differentiating  (23)  with  respect  to  X,  taking  appro¬ 
priate  care  of  the  dipole  singularity  in  the  volume  integral  at  the  pivotal  point  (X,  Y,Z).  The 
form  of  the  nonlinear  singular  integral  equation  for  U  depends  upon  the  definition  of  the  principal 
value  of  the  singular  volume  integral.  Following  Ogana  [12]  we  define  the  principal  value  of 
this  integral  by  surrounding  the  singularity  with  a  sphere  of  small  radius  and  take  the  limit  as 
the  radius  tends  to  zero.  It  follows  from  (23)  on  differentiation  with  respect  to  X,  that 

E(X)U(X,T)  =  £(X)f/2(X,7')/6+f/p(X,r) 


00 


f/2(Xi,7'-0)Ai(X,Xi)  +  -UHXuT-®)R2{X,Xi) 
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(2Aa) 
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(246) 

(24c) 

(24rf) 


! 

where 

Ki{X,Xi)  =  {(K-yi)*+(Z-2i)*-2(Jf-Jri)*}/|X-Xi|5. 

X2(x,Xi)  =  2A/(2r-2ri)/|x-Xi|3+|x-XitA:,(x.Xi), 

X*(X,Xi)  =  -[Af|X-Xi|-(2r-W/IX-Xi|» 
and 

l/p(X,r)  =  ^«>p(X.r).  (24e) 

In  (24a)  (/p(X,r)  is  uniquely  determined  from  (24e'),  (226)  and  (14)  once  the  shape  and  motion 
of  the  body  are  specified.  It  is  easily  seen  that  equation  (24a)  involves  only  U,  Ur  and  C/tt.  If 
we  specify  that  U  is  zero  everywhere  for  T  less  than  zero  and  start  the  body  moving  at  T  equals 
zero,  then  using  an  appropriate  numerical  approximation,  (24a)  can  be  used  to  step  forward 
in  time  evaluating  U  at  all  points  of  the  flow.  Equation  (23)  can  be  used  on  the  same  way  to 
evaluate  d>  and  <I>t,  while  the  pressure  coefficient  can  be  found  from  (22e). 

i 

) 

» 

5.  CONCLUSIONS 

We  have  obtained  an  integro-differential  equation  for  the  streamwise  velocity  component 
U,  which  asserts  that  the  disturbance  outside  the  wing  and  wake,  is  the  same  disturbance  that 
would  be  produced  by  a  fictitious  distribution  of  sources  on  the  wing  surface,  a  vortex  sheet 
on  the  Z  =  0  plane  and  a  volume  distribution  of  streamwise  directed  doublets.  It  should  be 
noted  that  if  the  flow  is  steady,  then  derivatives  with  respect  to  time  are  zero  and  (24a)  then 
reduces  to  the  formula  given  by  Ogana  [12]. 

Both  (23)  and  (24a)  are  suitable  for  fast  numerical  computations.  It  has  been  stated  that 
these  equations  may  be  solved  by  a  time-stepping  routine  which  evaluates  (/  or  <6  at  a  given  point 
of  the  flow  from  values  obtained  at  earlier  times.  It  should  be  noted  that  (23)  and  (24a)  may 
also  be  solved  by  assuming  that  the  motion  consists  of  a  steady  part  plus  a  small  harmonically 
oscillating  unsteady  part.  Thus  time-independent  integro-differential  equations  would  be  ob¬ 
tained  which  can  be  solved  by  an  iterative  technique. 
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APPENDIX 


lategral  EqaatlM  FomalatkM  for  Sqpenoaic  Mach  aanbers 

When  the  free  stream  Mach  number  is  supersonic  (M  >  I)  then  the  Green’s  function  for 
(4<i  and  Ab)  is  (see  Morino  [7]) 

C(x-xi,  t-h)  =  — (25a) 
where 

r,(x,xi)  =  (2») 

«*  =  A/2_l  = -/Ji!  (25c) 

and 

+  r.(x,xi)l.  (25<0 

a* 

In  (25^1,  256  and  25d)  it  has  been  assumed  that 

>  (x[(y-yi)^+(z-2iy^]K  (26) 

If  (26)  is  not  satisfied  then  G(x— xi,  /— ti)  =  0.  ’ 

An  integral  equation  for  the  reduced  velocity  potential  when  the  free  stream  Mach  number  is 
supersonic,  can  be  derived  by  following  a  similar  procedure  to  that  described  in  Section  2.  That 
is,  multiply  (8)  by  the  domain  function  £  (see  (6)),  substitute  for  G  from  (2Sa)  and  integrate  over 
both  space  and  time.  The  resulting  integral  equation  is  similar  to  (9)  (or  (10))  except  that  twice 
as  many  terms  appear  on  the  right  hand  side  due  to  the  two  time  delay  terms  in  (2Sa).  Note, 
as  in  the  subsonic  case,  the  contribution  from  shock  waves  will  not  explicitly  appear  in  the  integral 
equation  formulation  due  to  the  shock  capturing  nature  of  the  method. 

Applying  the  quasi  steady  approximation  as  derived  in  Section  3,  together  with  the  relevant 
boundary  conditions  (cf.  with  Section  3)  and  introducing  the  apposite  Prandtl-Galuert  variables 
(for  M>\)  (cf.  with  (12a  and  126)), 


X  =  X,  Y  =  ay,  Z  =  <xz  and  T  = 


M 


(27a) 


with 


O  =  M2(l +y*)^/a2. 


(276) 


we  find  that  the  perturbation  potential  satisfies  (cf.  with  (22a,  226  and  22r)  and  (136,  13c  and 
13c)), 


E(xyt>ix,T)  =  fl»p+(x,7^+<*>p-(x.r) 
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Cp  =  =  MH\  +y*)Cp/a2. 

(28c) 

0±  =  M(X-Xi)±R, 

(28rf) 

R  =  {(X-XiY-(  Y-  Ki)2-(Z-Zi)2}», 

(28e) 

dN-'  V  SX’dY’dZ)'  ’ 

(28/) 

A  =  \X-Xi\-{{  Y-  YtP  +  (Z-Zi)2}i 

(28g) 

and 


r»±  =  /(Xi,r-0*). 


(286) 


Here  as  in  Section  3,  ^  and  N  are  the  gradient  operator  and  unit  normal  to  the  surface  in  Prandtl- 
Glauert  space  and  H  is  the  Heaviside  step  function. 

Following  Section  4  it  is  consistent  to  assume  that  the  wing  is  a  thin  nearly  planar  body 
with  normal  direction  parallel  to  the  Z-axis.  Subsequently  we  find  that  equation  (28a)  reduces 
to  (cf.  with  (23)) 


E{Xyb(?L,T)  =  «I>p+(X.D+<I>p-(X,D 
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where  Rs  and  Rw  represent  the  projection  of  £b  and  Xw  onto  the  X—Y  plane,  and,  A  (  ) 
represents  the  value  on  the  upper  side  of  body  or  wake  minus  the  value  on  the  lower  side. 

The  appropriate  integral  relationship  for  the  streamwise  velocity  component  C/,can  now  be 
determined  by  differentiating  (29)  with  respect  to  X,  taking  due  care  of  the  dipole  singularity 
in  the  volume  integral  at  the  pivotal  point  (X,  Y^).  It  follows  from  (29)  then,  that  (cf.  with 
(24o-24e)) 

E(X)U(\,T)  =  £(X){/2(x,r)/3+ (/p+(X,r)+ (/p-(x,r) 
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where 


^i(X.Xi)  =  -2{2(jr-^i)a+(y-yi)s+(z-Zi)*}//i6, 

^8±(X.Xi)  =  -2Af(X-Xi)/R»±RRiiX,Xi), 

^8^X,Xi)  =  -[MR±(X-Xi)}ilR» 


(306) 
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and 


C/p±(X,r)  = 


(30c) 


As  in  equation  (24a),  (30a)  involves  only  U,  Ut  and  C/tt.  Then  using  an  appropriate  numerical 
approximation,  (30a)  can  be  used,  stepping  forward  in  time,  to  evaluate  V  at  all  points  of  the  flow. 
Equation  (29)  can  be  used  in  the  same  way  to  evaluate  d>  and  d>T,  while  the  pressure  coefficient 
can  be  found  from  (28c). 
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